Interaction corrections to tunneling conductance in ballistic superconductors 
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It is known that in tlie two-dimensional disordered superconductors electron-electron interactions in the 
Cooper channel lead to the negative logarithmic in temperature correction to the tunneling conductance, 
Sgoos -ln(^:^), above the critical temperature T^. Physically this result appears due to the density of 
states suppression by superconductive fluctuations near the Fermi level. It is interesting that the other correc- 
tion, which accounts for the Maki-Thompson-type interaction of fluctuations, is positive and exhibits strong 
power law, SgMT if^)^' which dominates the logarithmic term in the immediate vicinity of the critical 
temperature. An interplay between these two contributions determines the zero-bias anomaly in fluctuating su- 
perconductors. This paper is devoted to the fate of such interaction corrections in the ballistic superconductors. 
It turns out that ballistic dynamic fluctuations perturb single-particle density of states near the Fermi level at the 
energy scale e ~ s/TdT - Tc), which is different from e ~ T - T^, relevant in the diffusive case. As the con- 
sequence, fluctuation region becomes much broader. In this regime we confirm that correction to the tunneling 
conductance remains negative and logarithmic not too close to the critical temperature while in the immediate 
vicinity of the transition we find novel power law for the Maki-Thompson contribution, Sgmr °^ (t^t^)^''- We 
suggest that peculiar non-monotonous temperature dependence of the tunneling conductance may be probed via 
magneto-tunnel experiments. 

PACS numbers: 74.40.-n, 74.25.F-, 72.10.-d 



As it is well known,' the leading-order fluctuation correc- 
tions to the conductivity due to electron-electron interaction 
in the Cooper channel in the vicinity of the superconducting 
transition are given by the Aslamazov-Larkin^ (AL), Maki- 
Thompson^'* (MT), and density of states^ (DOS) contribu- 
tions. The first one has a simple physical meaning of the di- 
rect charge transfer mediated by fluctuation preformed Cooper 
pairs. The other two contributions have a purely quantum ori- 
gin. The Maki-Thompson process can be understood as the 
coherent Andreev reflection of electrons on the local fluctu- 
ations of the order parameter while density of states effects 
originate from the depletion of energy states near the Fermi 
level by superconductive fluctuations. The relative importance 
of these three contributions depends on whether a supercon- 
ductor is diffusive (Trei <K 1), ballistic (Trei » 1) or granular 
(5 <K r <K max{ETh, T}). Here rei is the elastic scattering 
time on impurities, 6 is the mean level spacing in the grain, 
r is the escape rate, Eth - Djfi is, Thouless energy for a 
grain with the typical size and D is the diffusion coeffi- 
cient. The Aslamazov-Larkin correction is essential in both 
pure and impure superconductors. The Maki-Thompson is 
important only in the disordered systems, since there is strong 
cancellation between MT and density of states effects in the 
ballistic regime.* Usually unimportant DOS contributions be- 
come crucial in the systems containing tunneling junctions^'^ 
or in granular superconductors.^'"^ Tunnel barriers or granu- 
larity require multiple electron scattering for AL and MT con- 
tributions to be important. As the result, the magnitude of 
these effects is suppressed either by an extra powers of tun- 
neling matrix element ~ \tpk^ (in the case of tunnel barriers^) 
or by the small ratio gilgt, <^ 1 between inter-grain qy ~ r/6 
and intra-grain ~ EjiJd conductances (in the case of gran- 
ular superconductors^'^). 

In the study^ of tunneling anomaly between diffusive thin- 
film superconductors separated by an insulating layer it was 



shown that there is one specific MT-type process that con- 
tributes significantly to the renormaUzation of the tunneling 
conductance. This process appears to the first order in tun- 
neling probability \tpk'^, like DOS contribution, however, to 
the second order in interaction, unlike DOS, thus contain- 
ing one extra power of small Ginzburg number, Gi < 1. 
The reason why these two contributions have to be accounted 
simultaneously is twofold. Unlike DOS correction dgoos, 
which leads to the suppression of the tunneling conductance 
above the critical temperature Tc, the novel MT-type contribu- 
tion 6gMT leads to its enhancement. Second, although being 
smaller by the inverse power of dimensionless conductance 
g - vD = kf(^\l2n » 1, where v - m/j: is the single-particle 
density of states in two-dimensions and 4i = vfTsi, this spe- 
cific MT contribution has much stronger power-law tempera- 
ture dependence, 6gMT Gi^ijSr)^ - opposed to the weak 
logarithmic in temperature correction coming from the den- 
sity of states, dgoos -Giln{j^). One should recall here 
that the parameter that controls perturbative expansion over 
superconductive fluctuations is set by the Ginzburg number, 
which is just inversely proportional to the dimensionless con- 
ductance, Gi oc l/g. So that, this is really a competition be- 
tween these two contributions that defines the nature of zero- 
bias anomaly in fluctuating superconductors. As a result, due 
to an opposite signs of dgoos and 6gMT terms the full conduc- 
tance correction dgoos + SgMT has non-monotonous tempera- 
ture dependence that even may change sign. Similar observa- 
tions emerge in the context of granular superconductors."^ 

The reason for such strong temperature dependence of 
SgMT was attributed in Refs. 7 and 10 to the importance of ver- 
tex renormalization by coherent impurity scattering (Cooper 
ladders). If so this would imply then that anomalous Maki- 
Thompson contribution is absent in ballistic tunnel junctions. 
Such conclusion is also appealing in the view of strong can- 
cellation of MT effects in ballistic regime of superconducting 
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FIG. 1: DOS and MT tunneling conductance correction diagrams. 



thin films. ^ However, as we show in this work, in contrast to 
the expectation MT interaction correction to conductance in 
balHstic tunnel junctions remains important. It also exhibits 
strong temperature dependence, similar to that in the diffusive 
regime, but with the fractional powers depending on dimen- 
sionality. 

In what follows, we carry out a microscopic calculation of 
interaction corrections to the tunneling conductance in ballis- 
tic superconductors Tt^\ » 1 with the help of standard tem- 
perature diagrammatic technique.'' Within this formalism the 
conductance 
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with d - 1, 2, 3 be- 



ing effective dimensionality of a superconductor {Id wire, 
2d film or 3(i bulk). This approximate form of the inter- 
action is obtained from the general expression' L '(^, w) = 
-v[ln ^ + ^(i + M) + ^\TT,,)q^ - where ^{Tt,\) = 

^[•^(5) + A^n^.'l^'il) - + 4;;f;-)]' ™der the assumption 
that characteristic energies of fluctuations are small as com- 
pared to the temperature, at ~ T - Tc T, and momenta are 
small as compared to the inverse thermal length {j = ff/T, 
q ~ y/l/BTGL ^ which allows to expand digamma func- 
tions tfr at small argument. 

Since matrix elements f^jt depend weakly on the momenta 
near the Fermi-surface one can substitute summation over p 
and k in Eq. (3) by the corresponding integration over the ener- 
gies Zpki- ■ ■) X d^pd^ki- ■ ■). Once these integrations 
are performed assuming ballistic limit max{e„, cjk} » t^} and 
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gr = -e— Im[n^(Q)]„^ 



(1) 



is determined by the retarded component of the polariza- 
tion operator n*(Q). Here e is the electron charge and V is 
the voltage applied across the junction. In the case of non- 
interacting electrons Matsubara version of the polarization op- 
erator is given by the simple loop diagram, which reads an- 
alytically as n(Q,„) = r2f„ I,pk\tpk\'^G(p,e„ + Q,„)G(fc, e„), 
where tpi; stands for the tunneling matrix element, e„ - 
2jiT{n+\ 12) and Q.„, - 2nTm are fermionic and bosonic Mat- 
subara frequencies, respectively, and 
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defines the single-particle Green's function. Under the as- 
sumption of momentum-independent tunneling amplitudes a 
simple calculation then gives for the bare value of the conduc- 
tance gr - ^e^v^\tpic\^. The first-order interaction correction 
is given by the diagram shown in Fig. la 

dYlDosi^m) ^2r-Y^Yj \tpkfG\p, en) 

eii<^k pkq 

G(k, e„ + Q.m)G(q - p, Uk - £i,)L(q, <^k) , (3) 



where 6(x) is the step function, one can complete summa- 
tion over the bosonic frequency cok in Eq. (3) by converting it 
into the contour integral and make an analytical continuation 
ie„ — > e -H iO. By combining the result for STIdos with Eq. (1) 
one obtains density of states type correction to the zero-bias 
conductance 
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(2) where we introduced Keldysh component of the interaction 
propagator LF{q, cS) - [Lf^iq, ai) - L'^iq, o))] coth ^ while the 
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retarded/advances components U^'^^Hq, <jS) are obtained from 
Eq. (4) by the replacement \(jL>k\ — > +iu). The most singular 
in r - Tc contribution to dgoos comes from the branch-cut 
of L^iq, to) where w) tanh term can be ignored. By 
taking U^iq, u) ~ -(32iT^ /7Tv)[(Bq^ + tgl)^ + w^]"' one can 
complete energy integration in Eq. (6) and find 
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which amounts an insertion of a single interaction line into 
one of the Green's function and coefficient of two accounts 
for two such possibilities. This is the density of states effect 
since upper part of the diagram is just a self-energy for the 
Gip, ek). The interaction propagator is defined as 



where is the second derivative of the digamma func- 

tion. The remaining q sum is dominated by the small mo- 
mentum transfer where argument of the digamma function 
can be taken as the constant i/'"(l/2) - -14^(3) since 
max{B^^, Vfq, tgl) ^ T . One obtains then as the result'^ 
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where dimensionless prefactors are Ci - 



-4-z and 

p-fS 



vS ^/BT^ 

oc with S being the cross-section area of 

the wire and h being the thickness of the film. The numeri- 
cal coefficients are ai » 2.17 and 02 ~ 0.17. For the bulk 
3ii junctions dgoos/gr ^ -C3 = is small and tempera- 
ture independent. Notice also that 'for 2t/-case C2 is linearly 
proportional to the Ginzburg number One sees from Eq. (8) 
that strong suppression in the density of states near the Fermi 
level translates only into moderate renormalization of conduc- 
tance dgoos ■ This observation brings us to the necessity to 
study contributions to conductance coming form the interact- 
ing fluctuations shown diagrammatically in Fig. lb. The rea- 
son for this is similar to that in the diffusive regime. First of 
all, this contribution is of the same order in tunneling ~ \tpk\^ 
as the density of states one. Fig. la. Second, although having 
an extra small prefactor, Cd, this contribution is positive, un- 
like dgoos 1 and has much stronger temperature dependence, 
which may dominate dgoos in the near vicinity of the critical 
temperature. The competition between these terms defines the 
nature of zero-bias anomaly in fluctuating regime of ballistic 
superconductors . 

The diagram in Fig. lb defines Maki-Thompson correction 
to the polarization operator, which reads explicitly as 

dUMji^ni)^ ^ G^iP^ en)G(q - p,Uk- e„) 

G^(k, e„ + £l„,)G(q' -k,(ij[- e„)L(q, LOk)L(q', w^) . (9) 

It is important to comment here that although this correction 
looks like second-order DOS, it in fact contains the mixture 
of advanced and retarded blocks of the Green's functions, 
which by its analytical structure is the same as in the Maki- 
Thompson diagram. This is precisely the reason why this 
term is strongly temperature dependent. DOS effects always 
involve Green's functions of the same causality and thus bring 
subleading temperature dependence. One calculates momen- 
tum integrals in Eq. (9) by the prescription defined in Eq. (5) 
and after the analytical continuation finds corresponding cor- 
rection to the conductance 



6g^ 



cosh^^iTco ^'^'^ 



gj 2n'^T 
Im[L«(^,w)]Im[L'*(^',w')] coth(f:)coth(ff ) 
{vf ■ q + CO - 2e+) (vf ■ q' - u)' + 2e_) 



(10) 



where e± - e ± /O stands as the reminder of analyticity. After 
the consecutive energy integrations this formula simplifies to 
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As compared to the corresponding result in the diffusive case^ 
the novel feature here is appearance of the Uf-q factors, which 
limits the phase space for the momentum transfer and in a 





FIG. 2: Aslamazov-Larkin interaction corrections to the tunneling 
conductance which appear in the higher order in transparency ~ \tpi;f 
then DOS and MT contributions shown in Fig. 1. 



way changes power-law behavior of the singular term in the 
conductance. The remaining momentum integration can be 
completed in the closed form and gives 



SgA 



= bdCliTcTGL)' 



(12) 



fort/= l,2,3withZ7i = 0.06, /72 = 1.4- 10-^ and ^73 =4-10-1 
As anticipated 5gMT is positive and has much stronger power- 
law temperature dependence then dgoos near Tc- We con- 
clude here that anomalous temperature dependence of 6gMT 
known from the impurity vertex renormalization in the diffu- 
sive case^ survives in the ballistic regime as well, however, 
with the fractional powers of T - Tc- 

There are two more diagrams in the second order in in- 
teraction that contribute to the conductance renormalization. 
These are shown in Fig. 2 and define Aslamazov-Larkin cor- 
rections. However, unlike DOS and MT terms in Fig. 1 these 
contributions appear only to the second order in the tunneling 
transparency and thus contain and extra smallness ~ \tpk\*- An 
estimate for these diagrams for two-dimensional case gives 
SgAi/gr 0^ Itpki^C^iTcTGhy^^, which has smaller amplitude 
and weaker temperature dependence then 6gMT in Eq. (12). 

There is simple physical picture which allows to under- 
stand these results at the qualitative level. The current in 
a tunnel junction is determined by the product of the den- 
sity of states convoluted with the difference of Fermi func- 
tion, which measure occupation of the given state, namely 
/(V) ~ / de[/f (e + eV) - /F(e)]v(e + eV)v(€). Within the 
linear response one can identify then from I(V) the zero-bias 
DOS: 6goos/gT ~ / ^ cosh-^ andMT: dgMr/gr ~ 

/ ^^^T^ cosh"^ (^)^ conductance corrections. Thus, estima- 
tion of the temperature dependence of dgoos and 6gMT re- 
quires knowledge of the detailed structure of the density of 
states above Tc- To this end, let us understand at which en- 
ergy window 6e superconductive fluctuations deplete single- 
particle energy states near the Fermi level and what is the 
depth of this suppression. The energy scale can be esti- 
mated knowing the time needed for the superconductive 
fluctuation to spread over the distance of coherence length 

^{T) - ^0 ^J^^■ In the disordered case is determined by 
the diffiisive motion of particles and gives for 6e via the un- 
certainty relation de ~ t^^ - D^^^(T) = t^\^ cc T - Tc- In the 

clean limit ballistic motion defines another scale'^ de ~ = 
Vf^^^iJ) cc y/TciT - Tc). The depth of the depletion region 
in DOS, 6v{e) - -^ImS*(e), follows from the self-energy of 
the electron Green's function S(e„) = T Yjpka>„ G^ip, e„)G(q - 
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p,(jiik - enW^(q,en,0Jk - ^n)L{q,u)k), where impurity vertex 
T{q,en,e^) = T^le{-e„eJI{Dq^ + |e„ - ej) is present only 
in the diffusive Umit. Having calculated S^(e) DOS renor- 
malization reads 5vD(E)ie) ^ (TcTgl)^ FD(B)('2.eTGL), where 
^D(B)(2eTGL) are energy depending scahng functions which 
are universal for the given dimensionality. For example, in 
the diffusive 2<i-case^ 
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, (13) 



which agrees with our explicit diagrammatic calculation 
[Eqs. (9)-(12)]. The reason why dgoos remains logarithmic in 
both cases is due to the conservation law 6vD(B)ie)d€ = 0. 
Indeed, when performing energy integration in Sgoos Iqt ~ 
J cosh"^ (^)^ necessarily accounts for the pole of 
the Fermi function which set the relevant energies to be of the 
order of T and not T - Tc- For e ~ T both scaling functions 
FD(B)i^) coincide to the leading singular order in T - T^. 



while in the balUstic regime^^ 
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where x - (jr' ll^{2>))TcTGh- The two basic properties of the 
scaling functions are Fuie ^ 0) ^ const while ^^(e 
0) ^ Ijx which is actually valid for any dimensionality and 
also F£)(B)(e)de = 0. The latter is the manifestation of 
conservation law for the total number of states. Knowing 
these facts one readily estimates (5vzj(0) oc {TcTc]^)~ and 

A—d 

6vb{0) oc (7'cTgl)~- The most singular contribution to the 
MT conductance renormalization comes from the energy re- 
gion of maximally depleted Sv{e) where interaction of super- 
conductive fluctuations is the strongest. The width of this re- 
gion is roughly 66 and, thus, interaction correction may be 
estimated as dgurlgr « Sv^{0)5e. For the diffusive case this 
gives 



-1 



(15) 



which reproduces results of Ref . 7 while in the ballistic case 



0^ (TcTgl) ^ (TcTgl) ' OC I - — — I , (16) 



gT 



T -Tr 



The possible way to probe these temperature anomalies in 
the conductance above Tc may be via magneto-tunneling. Let 
us recall that magnetic field H acts as an effective Cooper pair 
breaking factor that drives a superconductor away from the 
critical region. As the result, the relevant energies o) that de- 
termine conductance corrections in Eqs. (6) and (10) are set 
by the largest cutoflF between inverse Ginzburg-Landau time 
~ T - Tc and cyclotron frequency w// oc //, namely 
0) ~ max{a»H, Tq\]. So that by changing the field one effec- 
tively traces temperature dependence of 6g. Although an ex- 
plicit calculations of magneto-conductance in ballistic super- 
conductors is quite involved task we rely here on the plausible 
suggestion that is based on the results known for the diffusive 
case.^^ One may expect logarithmic in magnetic field depen- 
dence for the DOS correction dgoosiM)! Qt oc - \n{TcliOH), 
when (jJh ^ Tq^, and a power law of H for dguiiH). 
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